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Exact cosmological solutions for effective actions in D dimensions inspired by the tree-level superstring action are 
studied. For a certain range of free parameters existing in the model, nonsingular bouncing solutions are found. 
Among them, of particular interest can be open hyperbolic models, in which, without any fine tuning, the internal 
scale factor and the dilaton field (connected with string coupling in string theories) tend to constant values at late 
times. A cosmological singularity is avoided due to nonminimal dilaton-gravity coupling and, for D > 11, due to 
pure imaginary nature of the dilaton, which conforms to currently discussed unification models. The existence of 
such and similar solutions supports the opinion that the Universe had never undergone a stage driven by full-scale 
quantum gravity. 



1. Introduction 

It is widely recognized that the early Universe at sub-Planckian scales has been a scene for many strong effects 
which are nowadays either weak or even unobservable. Even remaining on a semiclassical level, one has to take 
into account the quantum properties of matter, the probable dynamical nature of extra dimensions predicted by 
modern unification theories (strings, supergravities, etc.), the existence of unusual kinds of matter such as strings 
and/or branes as well as modification of gravity at high energy densities and space-time curvatures. 

Many of these features show the existence of semiclassical and even classical (tree-level) mechanisms able to 
circumvent the well-known singularity theorems and to prevent the formation of a cosmological singularity, keeping 
the curvature on sub-Planckian scales. This means that very probably there has not been an epoch of cosmological 
evolution driven by full-scale quantum gravity. 

Thus, some nonsingular models have been built with the aid of gravitational Lagrangians nonlinear in curvature 
(e.g., |l], |j), whose origin may be explained by the quantum properties of matter fields. Supergravities and string 
theories, the candidate "theories of everything", also suggest new opportunities. Though, if the extra space-time 
dimensions, being an inevitable ingredient in such theories, are considered dynamically, the singularity problem 
becomes even more involved since, in addition to the usual cosmological scale factor, the extra dimensions can 
collapse or blow up, leading to a curvature singularity. 

In this paper we will deal with particular multidimensional cosmologies obtainable from an effective action 
which conforms to the low energy limit of some currently discussed unification theories — see Q and references 
therein. (For recent reviews on string cosmologies see 0.) This effective action admits as many as three natural 
mechanisms able to violate the usual energy conditions and hence potentially lead to nonsingular cosmological 
models. One mechanism is related to the nonminimal nature of the dilatonic scalar field and its interaction 
with antisymmetric form fields. Similar models in 4 dimensions have been discussed in Ref. 0. Second, the 
Brans-Dicke coupling constant to may have values leading to a "wrong" sign of the dilaton kinetic term in the 
Einstein- frame Lagrangian, which happens in dimensions larger than 11 Q (in other terms, the dilaton becomes 
pure imaginary). Third, in some field models of string origin, formulated in space-times with multiple time 
coordinates, antisymmetric forms can have Lagrangians with a "wrong" sign 0. 

The models to be discussed are quite simple but natural from the viewpoint of the underlying theories. In a 
sense, our "vacuum" approach is alternative to that of Ref. || where a hot brane gas is considered whereas global 
antisymmetric forms are ignored. We shall see that the first two mechanisms (but not the third one) really work 
and lead to globally regular cosmological solutions to the field equations. 

Our aim here is only to demonstrate the effect of these mechanisms, therefore we do not consider the totality 
of exact solutions that might be obtained for the action ([!]) but restrict our attention to the simplest case: a single 
dilatonic field with a Brans-Dicke type Lagrangian, a single antisymmetric form of axionic nature and a space-time 
with only two scale factors: external, a(t), and internal, b(t) . We assume the external 3-dimensional space to be 
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isotropic (spherical, flat or hyperbolic) and seek globally regular solutions such that, at late times, a(t) exhibits 
expansion while b(t) and the dilaton cp(t) tend to finite constant values. We will specify the requirement to the 
"input" parameters of the theory that lead to the existence of models with the desired properties. We shall see 
that they exist among closed (spherical) and flat models only for special values of the integration constants, in 
other words, require fine tuning. Unlike that, favourable hyperbolic models appear without fine tuning in a certain 
range of the input and integration constants and therefore seem to be much more realistic. 

There are many other arguments in favour of open cosmologies | fL0| . The observations are known to give the 
cosmological density factor £7 smaller or close to unity; meanwhile, the presently popular spatially flat cosmologies, 
most convenient for various calculations, require the precise equalty f2 = 1, actually a sort of fine tuning. It is 
much more probable that the real Universe at least slightly violates this special requirement. 



The model 



Consider the action of D -dimensional gravity interacting with a dilatonic scalar field $ and antisymmetric forms 
F s , F p , account the contributions from both the Neveu-Schwarz — Neveu-Schwarz (NS-NS) and Ramond-Ramond 
(RR) sectors: 

(d$) 2 



Si = I d D x^g\ * 



Tl-uj- 



$2 



where 1Z is the scalar curvature, 



E 



Vs 



Fi 



E 



Vr 



(1) 

1 , u> is a (Brans- 



9 = | detail, (d$) 2 = g MN d M ^d N ^>, M,N = 0. 
Dicke type) coupling constant, n s and n r are the ranks of antisymmetric forms belonging, respectively, to the 
NS-NS and RR sectors of the effective action; for each n-form, F 2 = F rii M 1 ...M n F^ Il — Mn ; the sign factors r] s = ±1 
and rj r — ±1 depend on a particular field model. 

The action % is written in the so-called Jordan conformal frame where the field $ is nonminimally coupled 
to gravity. This form is actually obtained in the weak field limit of many underlying theories as the framework 
describing the motion of fundamental objects, therefore we will interpret the metric gMN appearing in (|l|) as the 
physical metric. Thus, if the fundamental objects are strings, one has in any dimension to — — 1, while in cases 
where such objects are p-branes, one finds |3| 



(P-l)(p-l)-(p+l) 2 
(D-2)(p-l)-(p+l)3' 



(2) 



where p is the brane dimension and D is the space-time dimension. The NS-NS sector of string theory predicts a 
Kalb-Ramond type field with n s = 3; the type IIA superstring effective action contains RR terms with n r = 2, 4, 
while type IIB predicts n r — 3, 5. The action (|l|) may also represent the bosonic sectors of theories like 11- 
dimensional supergravity (where the dilaton is absent, and there is a 4- form gauge field), or 10-dimcnsional 
supergravity (there is a dilaton and a 3-form gauge field), or 12-dimensional "field theory of F-theory" admitting 
the bosonic sector of 11-dimensional supergravity as a truncation. The model Q contains a dilaton and two in- 
forms of ranks 4 and 5; it admits electric 2- and 3-branes and magnetic 5- and 6-branes. The "wrong" sign r) r = — 1 
is found in IIA* and IIB* supergravities, obtained with timelike T-duality from IIB and IIA theories, respectively, 
and also in the field limit of M* theory, appearing as a strong coupling limit of IIA* theory j|, [ll| . 
The standard transformation 

g M N = ^ D -^9MN (3) 
leads to a theory reformulated in the Einstein conformal frame, more convenient for solving the field equations: 



s E = [ d D x ^\n v ^) 2 -E^i e2XsVp s ~ E e2KvF " 

j ^ n s . n r . 



(4) 



where all quantities are written in terms of the Einstein- frame metric g MN ; <?e = | det g MN \ ; for the scalar field 
we have denoted 

D — 1 ( D - V 



D-2 



while the coupling constants X s and X r are 



n» - 1 



coxiD - 2) 
2n r - D 
2uj y {D - 2) 



; r\u = sign [ lj 

(NS-NS sector); 
(RR sector). 



D-2 



(•5) 



(0) 
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The sign factor ry w distinguishes "normal" theories = +1), such that the kinetic term of the p field in 
(|j) has the normal sign corresponding to positive energy, from anomalous theories where this sign is "wrong" 
(rju) = — !)■ It should be noted that many theories with D > 11 involve r) u = — 1. According to @, 



(D-2) 



1 



(£>-2)(p-l) 



(7) 



Evidently, under the condition (D — 2)(p — 1) > (p + l) 2 we have ^ = — 1. For p = 2, 5 this happens when 
D > 11 , and for p = 3, 4 when D > 10. 

The following table gives the values of to and rj^/tol for some particular space-time and brane dimensions. 
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oo 





12 


4 


-8/5 


-2 


10 


4 


2 


8/25 


12 


5 


-2 


-10/9 


10 


5 
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-2 


-9/8 


14 


2 


-4/3 


-4 


11 
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-5/2 


-18/25 


14 


6 


-16/11 


-132/49 


11 
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oo 





26 


3 


-17/16 


-48 


11 


6 


1/4 


36/49 


26 


4 


-50/47 


-1128/25 



Some comments are in order. First, the well-known result to = —1 for strings (p — 1) in any dimension is 
recovered. Second, one obtains to — oo for 2- and 5-branes in 11 dimensions, which conforms to the absence of 
a dilaton in 11D supergravity that predicts such branes. Third, in 12 dimensions one has rj^ = —1 for p < 7, 
and such a theory Ref. |j] does contain a pure imaginary dilaton: the F -forms of ranks 4 and 5 are coupled to 
a dilaton field ip with the coupling constants \\ = —1/10 and A2 = — Ai , respectively, while the product Xp is 
real. As is concluded in Ref. ||, for D > 11 "imaginary couplings are exactly what is needed in order to make a 
consistent truncation to the fields of type IIB supergravity possible" . In our (equivalent) formulation, p and A 
are real and the unusual nature of the coupling is reflected in the sign factor rj^ . 

Supersymmetric models with D = 14 are also discussed |lj, [TJ]], while D — 26 is the well-known dimension 
for bosonic strings. 



3. Solutions 

Let us specify the space-time structure and the Einstein-frame metric as 

M = R„ x M x Mi x • • • x M„, dimM 4 = d h 

n 

ds l = - e M«)^ + ^ e v W(is 2 



(8) 
(9) 



where it is a time coordinate ranging in l„cl and dsf are the u -independent metrics of the factor spaces Mj , 
assumed to be Ricci-flat for i — 1, . . . , n whereas ds^ in M describes a space of constant curvature K = 0, ±1 , 
corresponding to the three types of isotropic spaces; M is thus interpreted as an external (observed) factor space. 

There is a diversity of exact solutions for the action (|J) without L m in space-times like (^|), discussed, in 
particular, in Refs. [fj| [l(| (see also references therein). We will be only interested here in cosmological 
solutions for a very simple special case: a single antisymmetric form Fr^i from the NS-NS or RR sector, having 
a single (up to permutations) nontrivial component i*i,..d where the indices refer to Mo, and a single internal 
space Mi , so that in ([)]) i = 0, 1 , and ip = ip(u) . Then the field equations are easily integrated. 

Let u be a harmonic time coordinate for the metric (^), so that a = do/3 + di/3 1 . 

The F-form is magnetic- type; the Maxwell- like equations due to (|J) are satisfied trivially while the Bianchi 
identity dF = implies 

Fi...d„ = Qyfgo, Q = const, (10) 
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where go is the metric determinant corresponding to cLsq and Q is a charge, to be called the axionic charge since 
the only nonzero component of F can be represented in terms of a pseudoscalar axion field in do + 1 dimensions. 
The remaining unknowns are f3 , /3 1 and tp . 

In the Einstein equations TZ M — tjS^IZ — , written for the Einstein-frame metric (0), the stress-energy 



tensor has the form 



e 2a T^ = -ir 7j7 Q 2 e M ^ 1 + 2 ^diag(+l, [-l] do , [+l] d J - \^ diag(+l, [-l] do+dl ) (11) 

where the first place on the diagonal belongs to u and the symbol f]d means / repeated d times; rj F = ±1 is 
the sign factor of our _F-form, originating from rj s or r\ r in (|l|) or (|4). 

Due to the EMT property T" + = (where z belongs to Mo), the corresponding Einstein equation has 
the Liouville form a — /3° + Ko(do — l) 2 e 2a ~ 2/3 , whence 



do-1 



e^- a = S(-K ,k,u) d = { 



e ku , K = 0, keR; 

k -1 cosh/cu, K Q = 1, fc > 0; 

k- 1 smhku, K = -l, k> 0; (12) 

it, A'o = —1, fe = 0; 

k^ 1 sin/cu, T^Tq — ~~ 1, k < 0, 



where k is an integration constant and one more constant is suppressed by a proper choice of the origin of u . 
Eq. ( p"2] ) can be used to express (3° in terms of (3 1 . 

It is helpful to consider the remaining unknowns as a vector x A — (J3 1 , tp) in the 2-dimensional target space 
V with the metric 



<°">-('o* £)■ t«")-( 1/( J*' i). 



The equations of motion then take the form 



-rrfY^e 2 " (14) 



with the function y(u) = rfi/? 1 + Xip, representable as a scalar product of x A and the constant vector F in V: 

y(u) = Y A x A , Y A = (di, A), r A = (1/d, »fc,A). (16) 



Eq. (|15| ) is a first integral of (|14|) that follows from the (^) component of the Einstein equations. 
The simplest solution corresponds to Q = (scalar vacuum): 

/3 =c 1 u + c 1 , ip = c tp u + c tp , (17) 

where c 1 , c 1 , c v and c v are integration constants. Due to (|i~5|), the constants c 4 = (c 1 , c^,) are related by 

c A c A = dd^c 1 ) 2 + r^c 2 = —^—K. (18) 
v do - 1 

If Q 7^ 0, Eqs. ( |l4] ) combine to yield an easily solvable (Liouville) equation for y(u): 

jj + V F Q 2 Y 2 e 2 y = 0, Y 2 = Y A Y A = d x /d + Vu) \ 2 . (19) 

This is a special integrable case of the equations considered, e.g., in Refs. [jl4|, |ll| |l~6| |. 

Among the diverse solutions to (|l^) existing for different values of Y 2 , we will choose, for our purposes, the 
solution for Y 2 > 0. One of the reasons is that even for ?y w = — 1 one has Y 2 > for fields from the NS-NS sector 
in any dimension and for fields from the RR sector if D < 17. For Y 2 > 0, Eq. ( |l9| ) gives 

e -v(u) = \QjX_ cosh [/ l ( u + (20) 

where Y = lY 2 ] 1 / 2 , h > and u\ are integration constants. The unknowns x A are expressed in terms of y as 
follows: 

x A = ^y(u)+c A u + c A (21) 
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where the constants c A = (c 1 , c v and c A = (c 1 , c^) satisfy the orthogonality relations 

c A Y A = 0, 

Finally, the constraint (llq) leads to one more relation among the constants: 



c A Y A = 0. 



Y 2 



CAC 



do 



do-1 



K. 



(22) 



(23) 



4. Analysis of cosmological models 



4.1. Prelimaries 

In what follows, we put 
obtained in the above notations 



In what follows, we put do — 3, so that d\ = D — 4, and identify, term by term, the Jordan-frame metric ds 2 



exp 



2ip 



-du 2 



AS 2 (-K ,k,u) 



ds 2 



+ e^'dsf 



w 1 (D-2)J [2S(-K ,k,u) 
where the function S{., ., .) is defined in ([l^), with the familiar form of the metric 

-dt 2 + a 2 {t)dsl + b 2 (t)ds\, 



ds 2 = 



(24) 



(25) 



so that a(t) and b(t) are the external and internal scale factors and t is the cosmic time. 

To select nonsingular models, let us use the Kretschmann scalar K, = RainpqR MNP ® , which is in our case 
a sum (with positive coefficients) of squares of all Riemann tensor components Rmn ■ Thus as long as JC is 
finite, all algebraic curvature invariants of this metric are finite as well. For the metric (^) with do — 3 one has 
(the primes denote d/dt): 



K + a' z 



•di(di-l)-^ 



(26) 



By ( pq) , JC — > oo and hence the space-time is singular when a — * , a— >oo, fr^O or b — * oo at finite proper 
time £ . Accordingly, our interest will be in the asymptotic behaviour of the solutions at both ends of the range 
R« = («min, Umax) of the time coordinate u, defined as the range where both a 2 and b 2 in ( |25| ) are regular and 
positive. (Note that, as t — > ±oo, a singularity does not occur when b(t) — » 0, or a — > in case iCo = 0.) At any 
u £ M. u all the relevant functions are manifestly finite and analytical. The boundary values u m ax and u m in may 
be finite or infinite; a finite value of u max or u m i n coincides with a zero of the function (|l2|). 

Among regular solutions, of utmost interest are those in which a(t) grows while b(t) tends to a finite constant 
value as t — > oo . Any asymptotic may on equal grounds refer to the evolution beginning or end due to the 
time-reversal invariance of the field equations. We will for certainty speak of expansion or inflation, bearing in 
mind that the same asymptotic may mean contraction (deflation). 

Let us now enumerate the possible kinds of asymptotics. 



Type I: u — > ±oo , where 



dt 



3 (A-2fc)M 



3 AM 



,fl|«| 



(27) 



with k > and certain constants A and B depending on the solution parameters. A favourable asymptotic of 
a(t) takes place when A > 2k : 

(i) A>2k: t— >oo, a — t A K A ~ 2k ) (power-law inflation); 

(ii) A = 2k: t ~ |u| — >• oo, a ~ e kt (exponential inflation). 

A reformulation for k < is evident. The internal scale factor b(t) tends to a finite limit if B — 0, i.e., under a 
special condition on the model parameters (fine tuning). 

Type la: a modification of type I when k = , so that at u — * oo 



dt 2 



-3 n Au 



du 2 , 



.-1 „A 



e Au du 2 , 



(28) 
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If A > , we have, as desired, t — > oo and a — > oo ; the expansion may be called "slow inflation" since it is only 
slightly quicker than linear: the derivative da/dt ~ u, which behaves somewhat like Int. If A < 0, then a — > 
at finite i (singularity). As for b(t) , one may repeat what was said in case I. 

Type II: u — » 0, where the function ( |l2] ) tends to zero, so that S{— Kq, k,u) ~ u, while other quantities involved 
are finite. In this case 

dt 2 ~ l/u 3 , a 2 - l/u, b 2 -> fin. (29) 

According to (|2S|), i — > ±oo, a(i) ~ |i| (linear expansion or contraction), whereas both b(t) and </?(£) tend to 
finite limits since they do not depend on S{— Kq, k, u) . 

The dilaton (p in all cases behaves like ln&(i), but, in general, with another constant B in each particular 
solution. 

This exhausts the possible kinds of asymptotics for Y 2 > . Solutions with Y 2 < , which can emerge when 
rju = —1 and/or with ry F = — 1 , may have other asymptotics, but they are of lesser interest. 



4.2. Scalar-vacuum cosmologies 

The scalar- vacuum models (pi|), ( |l7|) depend on two input constants, D (or d\ — D — 4, or d — (D — 2)/2) and 
uj (or lji ) and three integration constants k, c 1 , c v related by (|l8|); two more constants, c 1 and , only shift 



the scales in Mi and along the (p axis and do not affect the qualitative behaviour of the models. 
Closed models, Kq = +1. In this case in (pi|) S = kj coshfcu, k > 0, hence the solution has two type I 



asymptotics at in ±oo , with k > and the following constants A = A± : 



d lC ! + 



(30) 



so that at least at one of the asymptotics A < whence a — » at finite £ , a singularity. It is also easily seen that 
if 6(t) const, it behaves as e Bu , B — const, and tends to zero at one of the limits u — > ±oo. 

Spatially flat models, iio = 0. One has simply 

a 2 {t)=e Au , dt~ e (A - 2k ^/ 2 du, (31) 

where A = —c v> /(duji) — die 1 — k, k 6 M, and again 6 2 (t) = c B ", i? = const. Thus each of the scale factors is 
either constant, or evolves between zero and infinity, and a — occurs at finite t. 

Hyperbolic models, Kq = — 1. If k > [note that, when rj u = 1, there is necessarily fc > due to (|l8|)l, one 
has in (24) S = k^ 1 sinhku. Hence the model evolves between a type I asymptotic at u — > oo, with A coinciding 
with A + in Eq. (|30|), and type II at u = 0. Since type II is regular, a necessary condition for having a nonsingular 
model is A > 2k. 

To find out if and when it happens for rj u = +1, it is convenient to introduce, instead of the two constants c 1 
and c v connected by (|8|), an "angle" 9 such that 



— c 1 = \ j, k cos 6*, — c v = y^fc sin#. (32) 

V 2ddi V 2 

The condition A > 2k will be realized for a certain choice of the integration constants if A + given by ([}0|) has, as 
a function of # , a maximum no smaller than 2k . An inspection shows that it happens if 

ujf < l/[d(6d - di)] = 1/[(D -1)(D- 2)]. (33) 

This is the only example of a nonsingular (bouncing) vacuum model with r] u = +1. 

In case k > 0, rjv = — 1, a choice of c v and c 1 subject to (|l^) such that A > 2fc is easily made for any wi . 
For riuj = — 1, /c = 0, the model evolves between type la and II asymptotics, where at the la end (u — > oo) 

A = -die 1 - c v /(duji), B = 2c 1 - c v /{du)i). (34) 

The necessary condition for regularity, A > , is satisfied for proper C\ and c v which can be chosen without 
problems. 

In case rj^ = — 1, k < 0, the function ( p"2| ) is simply Ifcl" 1 sin |fc|it, and the model has two type II asymptotics 
at adjacent zeros of S , say, u = and u = it /\k \ . This model is automatically nonsingular for any further choice 
of integration constants. 
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We conclude that among vacuum models only some hyperbolic ones can be nonsingular. For = +1 in such 
a case a(t) evolves from linear decrease to inflation, or from deflation to linear growth. Only in the latter case 
both b(t) and tp tend to finite limits as t — > oo without any fine tuning. 

For riu, = —I there is a model interpolating between two asymptotics of the latter kind. Thus, as t changes 
from — oo to +00 , a(t) bounces from linear decrease to linear increase (generically with a different slope) whereas 
b(i) and ip(t) smoothly change from one finite value to another. The latter model exists for generic values of the 
integration constants. 



4.3. Cosmologies with an axionic charge 

The solution contains, in addition to the input parameters D, lo and A, three independent essential integration 
constants: the "scale parameter" k, the charge Q and also h and c v connected by (^3|); the constant c 1 is 
excluded by the first relation ( p2|) 

die 1 + Xc v = (35) 

so that the quantity c a ca, appearing in (p3|), is expressed as c a ca = i]u(d/di)ctY 2 . The fourth constant, the 
"shift parameter" u\ , as well as c 1 and c^, connected by (22), are qualitatively inessential. 
Let us begin with "normal" models, rju> = +1- The solution ( f2l| ) has the form 

P\u) = ^y( u ) + c}u + c\ (36) 

<^( u ) = yiy( u ) + c v u + £ v - ( 37 ) 



with y(u) given by (|20|). The form of (|23|) implies k > and suggests a notation similar to (p2|), namely, 

3,,, „ L . , 2 d\~ 1/2 (3 



^kYcosB, c,p "{ 1+ d~) y7j ks ™6. (38) 

Let us now pass to the asymptotic description. 

Common asymptotic u —* 00 . This asymptotic is common to all K$ and belongs to type I with 

A= -k + hA 1 + Cip A 2 , B = hB 1 +c v B 2 (39) 



with the notations 



foAi + c v A 2 > 3fc, (41) 



The condition ^4 > 2/c , or 

required for nonsingular models, may be realized for small Using the representation (|38|), one can find, 

precisely as in the scalar- vacuum case, a condition under which j4 max j the maximum value of A as a function of 
8, satisfies ^4 max > 2fc. Curiously, the coupling A drops away from the calculation, and the resulting condition 
coincides with (|3^). 

Meanwhile, B can have any sign, therefore the asymptotic of b[t) is uncertain; however, by choosing the ratio 
c v /h (that is, by fine-tuning 0) one can achieve B = 0, so that b(t) has a finite limit. 
The solution behaviour at the other end of the range R„ depends on Kq . 

Closed models. For Kq = +1, u — * 00 the asymptotic behaviour is 

a 2 (t) ~ e A 'l"l, b 2 {t)~ e B ' |u| , e (A '- 2fe)|u|/2 ^ (42) 



where A' and coincide with A and B given by (|39j) with the replacement c v — > — c v (or 6* — > — ^). Therefore 
the description is the same up to the replacement t — > — t. When simultaneously A > 2k and A' > 2k, one 
has a deflation — > inflation transition for the scale factor a(t) , with generically different powers A/ (A — 2k) and 
A' j (A' — 2k) at the contraction and expansion phases, and a regular bounce between them; at one end the evolution 
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may be exponential. But, even if b(t) — > fin as t — > oo , the evolution of 6(i) begins with & = or & = oo unless 
the model parameters are further fine-tuned. 

Spatially flat models. The case Kq = differs from Kq = +1 in that k can have either sign, and when 
comparing the asymptotics u — ► oo and u — > — oo, one has to replace fc — > — fc. Hence, if we deal with a 
solution such that the Universe evolves in a power-law inflationary regime in the asymptotic u — > oo , then, in 
the asymptotic u — ► — oo, we find a behaviour like a(t) ~ |i| s , s < I. So in the remote past there might be 
a subluminal contraction while in the remote future the model inflates. A necessary condition for this type of 
behaviour is again ( |4l| ) (with k replaced by \k\), leading to the requirement (|33|). 
Concerning b(t) , one only can repeat what was said about K = +1 . 

Hyperbolic models. The second asymptotic, u — > 0, is type II. Thus, provided A > 2k (sec (^)), there occurs 
linear contraction of the Universe in the remote past and inflation in the remote future, or, in a time-reversed 
model, deflation in the remote past and linear expansion in the future. The latter opportunity seems especially 
attractive since both ip (hence the string coupling) and b(t) (hence the effective gravitational constant) tend to 
finite limits automatically, without any need for fine tuning. 

We see that even models where the kinetic terms of tp and Fmnp have both normal signs (rj^ — rj F = 1), 
predict some nonsingular bouncing cosmologies. 

"Anomalous" models with rj u = —1, like their vacuum counterparts, bear some new features compared to 
r) u = +1, connected with the relation (|2^). Namely, Eqs. (p0|), ([56]), (|37|) are still valid, but now we can have 
c a ca < 0. Therefore, first, the representation ( |38| ) is no more valid, thus cancelling the restriction (|33|). Second, 
for Kq — — 1 , one can have now solutions with k < 0, with a different analytical form and different behaviour. 

As a result of the first of these circumstances, one obtains a bouncing behaviour of a(t) in closed (Kq = 1) 
and flat (Kq = 0) models, having two type I asymptotics as described above, by simply choosing the integration 
constants h and ipi from a proper range, without further restrictions on the input parameters (provided they lead 
to rju — — 1 ). However, fine tuning is still necessary for obtaining finite limits of b(t) and (pit) at late times. 

The appearing hyperbolic models with k = interpolate between type II (u — > 0) and type la (u — > oo) 
asymptotics (linear contraction — > slow inflation or slow deflation — > linear expansion), under the condition 

A = hAx + c v A 2 > 0, (43) 

easily satisfied by choosing proper h and c v . 

There also appear bouncing models with k < which qualitatively behave quite similarly to their vacuum 
counterparts, interpolating between two type II asymptotics. 



5. Concluding remarks 

The bouncing mechanism discussed here works in a certain range of the free parameter of the model, the Brans- 
Dicke type constant u>. Part of this favourable range, see Eq. (|3^), corresponds to "normal" dilatonic fields with 
positive energy in the Einstein frame, and such bouncing models exist for all D > 4. If we, however, ascribe the 
origin of lo to fundamental p-branes in the spirit of Ref. ||, it turns out that only "anomalous" theories with 
rjuj = —I lead to bouncing models. [Indeed, in case rj w = 1 Eq. ([?]) gives 2 < D — 2, contrary to (J33|) - ] Moreover, 
such models only exist for D > 11 , see the table in Sec. 2 and the comments after it. 

The existence of a type II asymptotic of open models, making it possible to avoid fine tuning in getting a desired 
large t asymptotic, is quite a general phenomenon for cosmologies with the structure (H). It actually follows from 
a "stiff" character of the stress-energy tensor (pressure in the external space is equal to energy density) that leads 
to Eq. (|l2] ) and is common to the dilaton and the axion in the Einstein frame. 

The present nonsingular vacuum and axionic models evidently cannot pretend to describe the full-time evolution 
of the Universe but rather the epoch of maximum contraction under the assumption that this stage is dominated by 
scalar-vacuum and axionic effects. Unlike other scenarios with a bounce at sub-Planckian scales, such as the brane 
gas |9) and Pre-Big Bang || |T^| scenarios, no higher-order curvature terms are needed here to complement 
the model. Furthermore, our models, with b(t) and tp(t) tending to finite constant values, become effectively 
4-dimensional at late times, with constant values of the effective gravitational constant and string tension, so 
the further evolution can be studied using conventional methods, in particular, an inflationary period can follow. 
Though, in our view, inflation is mostly needed for solving the problems of Big Bang cosmology emerging due to 
the existence of multiple causally disconnected regions, whereas bouncing cosmologies do not create such problems. 

In any bouncing model a separate problem is the origin of its initial state. We would note here that in some 
cases, as i -> —00, the dilaton field $ = e - ^/" 1 tends to zero, so that the string coupling parameter g s = 1/$ 
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diverges. This problem may probably be coped with due to the superstring dualities, mapping a strong coupling 
regime to a weak coupling regime. Such a difficulty is, however, absent in the above models interpolating between 
two type II asymptotics. 
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